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le the earlier Parts 1 and 2 of the series (Bibliography, Nos. 4 and 5) a 
theory of the directly heated thermistor has been given. In the present Part 4 
this theory is extended to indirectly heated thermistors in order to make 
possible the handling of circuit problems where the dynamic properties of the 
indirectly heated thermistor are decisive. 

In particular a small signal circuit for the indirectly heated thermistor is 
derived. The parameters of this circuit follow from an analytic expression 
of the resistance-power characteristic, which in its turn can be deduced from 
a few simple measurements. 

The validity of the theory has been tested by experiments. The agreement 
is satisfactory. 
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List of Symbols 


¢ = time. 

Q = thermistor temperature. 
QO, = ambient temperature. 
U 

I 


= total thermistor voltage. 
= total thermistor current. 
U, = total heater voltage. 


I, = total heater current. 

uw = superimposed variable thermistor voltage. 
i = superimposed variable thermistor current. 
u, = superimposed variable heater voltage. 

i, == superimposed variable heater current. 


E, = emf on heater side. 
E, =emf on thermistor side. 


N = UI = power supplied to thermistor. 

N, = U;, I, = power supplied to heater. 

P = dissipated heat power. 

W =thermal energy of thermistor body. 

R =U /I = static resistance of thermistor. 

r = (dU/d1)y, — const, = dynamic resistance of thermistor. 
R, = U,/I, = heater resistance. 


R, = series resistance on heater side. 
R, = series resistance on thermistor side. 


Z = thermistor impedance. 

= — (€AN/dN,)p const. = heater efficiency. 

F N/R-dR/dP=d ic facto 
| = — = dynamic factor. 

adN/N Np, = const. 

m = 1/R-dR/dP = power coefficient of thermistor resistance. 

Rk, = constant in the resistance equation (ohms). 

B =constant in the resistance equation (°K). 

C = power sensitivity (°K/mW). 

T = thermal time constant. 


R, — Ry 


t =f7 ( + F = effective time constant of thermistor. 


Index 0 indicates value in operating point. 
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CHAPTER 1 


Introduction 


An indirectly heated thermistor consists of a small bead of ther- 
mistor material (“the thermistor”) with a large, negative temperature 
coefficient, surrounded by a heater with a very small temperature 
coefficient. The heater is electrically insulated from the thermistor 
by a thin layer of insulating cement. The device is enclosed in a 
small, evacuated or gas-filled glass bulb, see Figs. I to 3. 

The resistance of the thermistor is determined by its temperature, 
which can be changed by supplying current to the heater, to the 
thermistor or to both. Of course, the thermistor temperature also 
depends on the ambient temperature. On account of the small di- 
mensions of the thermistor, already a small amount of supplied 


Fig. 1. Indirectly heated thermistor in evacuated glass bulb, Stantel type B 2552/60, 
approximately full size (manufacturer: Standard Telephones and Cables Ltd). 
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Fig. 2. Schematic picture of indirectly heated thermistor. The numbers refer to Fig. 3. 


electric power gives a considerable temperature increase. Thus, a 
supplied power of only 3 mW increases the temperature of the ther- 
mistor in Fig. 1 with about 20°C and decreases its resistance to 
less than half the room temperature value. 

This thermistor type has a maximum allowable dissipated power 
of 60 mW, corresponding to a resistance less than 1/1000 of the room 
temperature resistance. The temperature coefficient of the heater, 
on the other hand, is quite low, so that for normal operating tempera- 
tures the heater resistance varies only a few percent. 

Great sensitivity, small dimensions and considerable sturdiness 
have made the indirectly heated thermistor widely employed in many 
fields [1, 2, 3,6]. For a long time it has been used in applications 
where its properties as an electric circuit element are of decisive 
importance. Nevertheless, a practically useful theory of circuits 
containing indirectly heated thermistors has not yet been published. 

In the earlier published Parts 1 to 3 of this study of thermistor 
circuits [4, 5, 7], the theory of the directly heated thermistor has 
been treated. The only fundamental difference between the directly 
and the indirectly heated thermistor is that the temperature of the 
indirectly heated thermistor can also be influenced through a heater. 
Therefore, the theoretical treatment of the two types has the same 
main features. The directly heated thermistor is equivalent to an 
indirectly heated thermistor, supplied with constant heater power. 
Hence, the results of this paper can also be immediately applied to a 
directly heated thermistor. 

In Chapter 2 the static properties of the indirectly heated thermistor 
are treated. It is shown that a few measurements are sufficient for 
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Fig. 3. Symbol for indirectly heated thermistor. 


determining the data necessary for the solution of a static problem. 
In Chapter 3 a derivation is given of the so-called thermistor equa- 
tion, which governs the dynamic behaviour of the thermistor. An 
example shows the application of the thermistor equation to the 
solution of a dynamic problem. An equivalent small signal circuit 
of the indirectly heated thermistor is derived in Chapter 4. Some 
applications are given. The validity of the theoretical developments 
has been tested by an extended series of measurements, treated in 


Chapter 5. 


The work described in this paper was made possible by grants 
from the Swedish Technical Research Council and from the Chalmers 
Research Fund. The author is indebted to Professor S. EKELGOF, 
Head of the Institute for Theoretical Electricity and Electrical 
Measurements, for having directed his interest towards the present 
problem and for valuable advice in the course of the work. He also 
wants to thank Tekn. lic. R. Davirpson for stimulating discussions. 
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CHAPTER 2 


Static Conditions 


a. Static characteristics 


Throughout this paper the ambient temperature and other external 
conditions, which influence the thermistor temperature, are supposed 
to be constant. 

The properties of the indirectly heated thermistor under static 
conditions are then completely defined by its static characteristics. 

The static U-I-characteristic (Fig. 4) gives the relation between 
thermistor voltage U and thermistor current J with the heater current 
I, or the heater power N), as a parameter. 
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Fig. 4. Static U-I-characteristics with the heater power Ny), as a parameter. 
Thermistor Stantel type B 2552/60 at O, = 25° C. 


a) Nz, = 0 d) N;, = 30 mW 
b) N, = 10 mW e) N;, = 40 mW 
c) Nz, = 20 mW 
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Fig. 5. Static R-N-characteristics of the thermistor of Fig. 4. 
a) Nz = 0 d) = 30mW 
b) Nz, = 10 mW e) = 40 mW 
c) Nz, = 20 mW 


The static R-N-characteristic (Fig. 5) gives the relation between 
the thermistor resistance R and the supplied thermistor power NV 
with J, or N; as a parameter. Evidently, each of these characteristics 
can be calculated from the other by the relations 


Ra Ol (1,2) 
or, inversely 
U = (NR); I= (N/R) (3,4) 


If we have an electric network containing a thermistor, we can 
use the static characteristics for calculating graphically in the usual 
way its currents and voltages. 

If high accuracy is desired, we must sometimes pay regard to the 
fact that in spite of its small temperature coefficient the heater 
resistance R, increases noticeably with increasing thermistor tempera- 
ture. It follows that the heater power NV, = R, J} takes on different 
values in the different points of a static U-J-characteristic even if 
the heater current J, is kept constant. Hence the static characteristics 
with the parameter NV, differ somewhat from those with the parametcr 
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I,. Generally, the variation in R, is so small that, for all practical 
purposes, it can be neglected (for the thermistor in Figs. 4 and 5 
the variation is less than + 1 percent). 

The procedure for determining the static characteristics is obvious. 
If we want the heater current as the parameter, we send a constant 
current J, through the heater and a constant current J through the 
thermistor. We then measure the voltage U when it has reached its 
steady state end value. The measurement is repeated for a second 
value of the thermistor current J and so on. If we want the heater 
power as the parameter, we must also measure the heater voltage 
U,, and adjust at each measurement the heater current J; so that 
N, = U, 1, = const. From the measured values of U and J we can 
directly plot the U-I-characteristic and, by the aid of (1) and (2), 
also the R-N-characteristic. The case N, = const. also gives us 
the variation of R, = U;,/I). 


b. Simplifying assumptions. Heater efficiency 


Unfortunately the task of determining a series of static charac- 
teristics like those in Figs. 4 and 5 is quite time-consuming. Good 
accuracy requires many measuring points and in going from measuring 
point to measuring point one has, on account of the thermal inertia 
of the thermistor, to wait a considerable time for the steady state. 
Every possibility to reduce the number of measuring points is therefore 
of value. In this respect we find great help in the following schematic 
picture of the thermistor, similar to the picture of the directly heated 
thermistor which has served as a basis for the earlier papers of this 
series: 


All parts of the thermistor body proper, i. e., the thermistor bead with 
its insulating cover and that part of the heater which is in immediate 
contact with the cover, have one and the same temperature O. The remain- 
ing parts of the thermistor, i. e., connecting wires, glass bulb etc. have 
the temperature O, of the surrounding medium. 


From these assumptions follows that, at a given temperature 0,, 
the power P, dissipated as heat from the thermistor body to the 
surroundings, is uniquely determined by the temperature 0. The 
resistance R of the thermistor bead being also uniquely determined 
by 90, we must have a one-to-one relation between R and P. This 
relation is the resistance-power-characteristic of the thermistor. 
Hence, to each value of the resistance R always corresponds one and 
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the same value of the dissipated heat power P. It follows that P 
is constant on any straight line through the origin of the U-J-diagram 
(resistance line, see Fig. 4); it is also constant on a horizontal line in 
the R-N-diagram (Fig. 5). It should be observed that the R-N- 
characteristics are static while the R-P-characteristics are also valid 
under transient conditions. 

In the stationary state the thermistor body is in thermal equilibrium. 
This implies that the dissipated heat power P equals the total electric 
power supplied to the thermistor. This total power consists of two 
parts. From the thermistor circuit the electric power N is brought 
to the thermistor body. But the thermistor body also becomes heated 
from the heater circuit. Of the entire heater power N,, however, a 
certain amount is lost in those parts of the heater which have no 
direct contact with the thermistor body. Only a certain part JN), 
of N,, is brought to the thermistor body. The factor » < 1 is of para- 
mount importance in our theory. We will call 4 the heater efficiency. 

Thus, the steady state power balance equation of the indirectly 
heated thermistor can be written 


P=WN-+ (5) 


If the approximations, on which we have founded our discussion, 
are physically reasonable, we should expect a constant value of the 
heater efficiency. The practical importance of this conclusion calls 
for a thorough investigation of its validity. This investigation, 
reported in Chapter 5, shows that, for the thermistor types investi- 
gated, the heater efficiency is practically constant so that it does not 
depend on the operating conditions of the thermistor. We hence put 


74 = constant (6) 


We can compute 7 from two pairs of power values, (N,, Nj) and 
(N., Ny2), corresponding to one and the same value of P (and of 
R). Eq. (5) gives us 


Ni + = No + (7) 


from which 


(8) 


! 
Niaz Nin 


0 20 40 60 


P mW 


105 


Fig. 6. Approximately linear relation between R;, and P for the thermistor of Fig. 4. 


Going to the limit of small changes in N and N,, we find from 


eq. (8) 
aN 


aN. (9) 
dN), 

This expression can be used as a general definition of the heater 
efficiency, valid also if 7 varies. In words: 


The heater efficiency is equal to the ratio between small changes in 
thermistor power and heater power chosen so that the thermistor resistance 
remains unaltered. 


On account of the fact that we have found 7; = constant, we need 
not determine all the U-J- or R-N-characteristics through direct 
measurements. It is sufficient to determine 7 and one single charac- 
teristic, the R-P-characteristic of the thermistor. Eq. (5) shows 
that P = N for N,=0. Hence, the important R-P-characteristic 
is identical with that static R-N-characteristic which corresponds to 
N, = 9. 

We have mentioned earlier that if we want to make very accurate 
calculations, we must also know how the heater resistance R,, varies 
with the operating conditions of the thermistor. Our schematic 
picture of the thermistor implies that that part of the heater which is 
in direct thermal contact with the thermistor bead has the temperature 
O of the bead, while the rest of the heater has the temperature of the 
surroundings. This means that we should have a one-to-one relation 
between R, and @ and thus also between R, and P. The investiga- 
tion in Chapter 5 shows that such a relation between R, and P can 
be assumed even if we aim at high accuracy (see Fig. 6). The relation 
being linear, its construction only requires that we know two points 
on the #,-P-characteristic. 
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c. The resistance equation 


In Part 2 of this study it has been shown that the R-P-characteristic 
of a directly heated thermistor can be represented analytically with 
good accuracy by the resistance equation 


B 
B. +6 (10) 


R,, (unit: ohm) and B (unit: °K) are characteristic constants of the 
thermistor in question, while the power sensitivity C (unit: °K/mW) 
also depends on the properties of the surrounding medium. Now, we 
know that the R-P-characteristic of the indirectly heated thermistor 
is obtained as the relation between resistance and thermistor power 
with the heater power = 0. An indirectly heated thermistor with- 
out heater power being essentially equivalent to a directly heated 
thermistor, it can be surmised that eq. (10) should be valid also for 
the indirectly heated thermistor. For the thermistors investigated 
in Chapter 5 this proves to be the case. 

By introducing (5) in (10) we obtain the expression of ‘the static 


R-N-characteristics: 
B 
R=R (11) 


Thus, the heater power N, has the same influence on the static 
characteristics as an increase of the ambient temperature by the 
amount C7N,,. 

It has also been shown in Part 2 how the three constants R, , B, 
and C can be calculated. The simplest method is the following. 
Starting with three points (R,, P,), (R,, P,), and (R,, P,) on the 
R-P-characteristic, we obtain from eq. (10) the two relations 


BC (P, — P 
In R, — In R, = ee (12 b) 


(CP, + po) (CP; + 


which give us 


(P, — P;) In R, + (P, — P,) n R, + (P, — P,) Ink, 13) 
= — Op (p,—P,) nR, + P,(P,—P,)InR, + P,(P,—P,) | 
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Introducing this value of C into (12a) or (12b), we obtain B. 
R,, then follows from (10) if we make use of one of the three known 
points, for instance (R,, P). 


From what has been said, we can infer that the problem of deter- 
mining the characteristics of the indirectly heated thermistor can be 
considerably simplified. In fact, the following measurements are 
sufficient: 


1) The R-P-characteristic is determined as the relation between 
resistance and dissipated thermistor power with the heater power = 0. 
If we use the resistance equation (10), we only need three points on 
the R-P-characteristic in order to know it completely. 


2) The heater efficiency is determined from two pairs of corres- 
ponding values of thermistor power and heater power, belonging to 
the same resistance value. 


3) The linear relation between FR, and P is determined by measuring 
R, for two different values of P. 


These measurements give us all the information about the in- 
directly heated thermistor which we need for the treatment of static 
problems. From the data they yield, we can compute in a simple 
manner the R-N-characteristics as well as the U-I-characteristics 
with N, as a parameter. R, being known, we can also obtain the 
characteristics with the parameter /,. 


d. Dynamic resistance, dynamic factor 


The small signal equivalent circuit, deduced in Chapter 4, makes 
use of the dynamic resistance r and the dynamic factor F of the ther- 
mistor. Both quantities are a measure of the derivative of the static 
characteristic in a given operating point. The following expressions 
are generalizations of the corresponding definitions in Part 1: 


(ar) 4) 

dl = constant 

N (=) ae 
dN |N = constant R dN Nj, = constant 
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Thus, the dynamic resistance is equal to the derivative of the 
static U-I-characteristic with NV, as a parameter, while the dynamic 
factor is proportional to the derivative of the static R-N-characteristic 
with N, as a parameter. 

Now, the two characteristics depend on each other. In fact, we 
have the following relation between r and F, already derived in Part 1: 


This can also be written 
2 1—F 
r= (17) 


In the maximum point of the static U-J-characteristic we have 
r= 0. It follows from (16) that in the maximum point F = 1; to 
the left of the maximum point we have r > 0, F < 1; to the right of 
it we have r < 0, F >1. In the origin we have F = 0 and r= R. 

In eq. (15) we have N, = const. Using (5) and (6) we can therefore 
write 


N 
(18) 


or 
F=— Nm (19) 


where we have introduced the power coefficient m of the resistance: 
C= (20) 


The power coefficient is equal to the ratio between a certain relative 
change of the resistance and the corresponding absolute change of 
the dissipated heat power. 

By the aid of eq. (10) we can obtain an analytical expression of 
the dynamic factor F, which obviates the need for graphic determina- 
tions of derivatives. Taking the logarithm of (10), and using (18), 
we immediately obtain 


BCN 
~ (OP + Oy) 


(21) 
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Fig. 7. Dynamic factor F as a function of N. Same thermistor as in Fig. 4. 
a) Nz = 0 d) N;, = 30 mW 
b) N, = 10 mW e) N;, = 40 mW 
c) Nz, = 20 mW 


or, after introducing (5): 


BCN 
= TON + On, + 0) (22) 


It follows that the dynamic factor at constant heater power has a 
maximum value 


B 
occurring for 
N = 9N, + (24) 


We see from (23) that the largest possible value of F is B/40), 
occurring when NV, = 0 and N = @,/C. For the thermistors investi- 
gated numerically in this paper, the constant B is of the order of 
magnitude 3 600° K (high resistance thermistors) or 2 400° K (low 
resistance thermistors). Hence, at room temperature (300° K) the 
maximum values of F are of the order of magnitude 3 and 2. For 
the thermistor of Fig. 4, Fig. 7 gives the relation between F and N 
with N,, as a parameter. 

When F has been determined from (22), r follows at once from (17). 
Fig. 8 gives the relation r vs. J with N;, as a parameter. 
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Fig. 8. Dynamic resistance r as a function of J. Same thermistor as in Fig. 4. 


a) N, = 0 d) N, = 30 mW 
b) N, = 10 mW e) N;, = 40 mW 
c) Nz, = 20 mW 


A comparison with the earlier papers of this series shows that our 
present formulae are immediately applicable to a directly heated 
thermistor, if we put V, = 0. 

Inversely, on account of the equivalence between an increase of 
the heater power and an increase of the ambient temperature, the 
earlier results are directly applicable to the indirectly heated ther- 
mistor. We only have to replace everywhere the ambient temperature 
0, with the “equivalent ambient temperature” 0, + CyN,. 

N. B. For a directly heated thermistor in the static case it is not 
necessary to distinguish between P and N. Our earlier formulae 
therefore sometimes have P where the present paper has N. 
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CHAPTER 3 


Dynamic Conditions 


a. The thermistor equation 


In the preceding chapter we have given the power balance equation 
of the thermistor under static conditions: 


N+7N,=P (5) 


Under dynamic conditions part of the supplied electric power 
serves to increase the thermal energy W of the thermistor body. There- 
fore, in this case the power balance equation reads 


N +N, =P +—— (25) 


In words: supplied electric power = dissipated heat power + in- 
crease per unit time of the thermal energy. 

As in the theory of the directly heated thermistor we introduce 
the thermal time constant T of the thermistor, 


26 
Then (25) can be written 
dP 
N+ 9N,=P+T— (27) 


This is the thermistor equation, generalized to the case of the in- 
directly heated thermistor. As shown in Parts 1 and 2, 7' varies 
slightly with the temperature 0 and thus also with P. Hence, in prin- 
ciple, 7’ in eq. (27) is not a constant but a function of P, 7’ = T (P). 
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0 20 40 60 


Fig. 9. Thermal time constant 7 as a function of P for the thermistor of Fig. 4. 


Fig. 9 gives the relation 7’ = 7(P) for the thermistor of Fig. 4. 
We see that 7' decreases with increasing P. The variations are small, 
however, so that it is often possible to put 7’ = const. 

Furthermore, the approximation involved in our assumption of a 
homogeneous temperature of the thermistor body makes it reasonable 
to neglect variations in 7’. In fact, especially under transient condi- 
tions only moderate accuracy can be expected from the thermistor 
equation. The heat transport within the thermistor body is not 
instantaneous, and hence a sudden increase of the heater power 
cannot immediately influence the temperature of the thermistor bead. 
For the same reason, a decrease of the thermistor power gives at the 
beginning of the transient period a higher temperature to the 
thermistor bead than the one deduced from our simplified theory. 

For the directly heated thermistor it has been shown in Parts 2 
and 3 that the uneven temperature distribution occasionally gives 
large deviations between theory and actual phenomena. If the 
power changes are not too large, however, the accuracy of the ther- 
mistor equation is sufficient. Hence, it does not seem necessary to 
develop a theory, which takes into account the uneven temperature 
distribution. 


b. Determination of the thermal time constant 


Dynamic problems can be solved by a direct use of the thermistor 
equation (27). If the thermistor current is so small that we can put 
N = 0, the calculations become especially simple. In such cases the 
thermistor equation can be written 


dP 


T sec. 
i 
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Fig. 10. Circuit for measuring the thermal time constant of an indirectly heated 
thermistor. 


With N,(t) known, we compute P(t) by direct integration and 
then obtain R(t) from the R-P-characteristic. Let us suppose that 
N, makes a sudden step change from the value N,, to the value 
Ni2- Then eq. (28) gives us 


P= + 9 (Naz — Na) (1 — (29) 


Thus, the dissipated heat power varies exponentially with the time 
constant 7 from the value P, = 7N,, to the value P, = 7N,.. This 
result gives us a simple experimental method for determining the 
thermal time constant, based on the circuit of Fig. 10. 

In this circuit the heater is connected in series with a resistance 
R,, which can be short-circuited by a switch S. The resistance R, 
should be much smaller than R,, so that the closing of S gives a very 
small change of N,. Then the thermistor resistance R also. changes 
very little and we can assume a linear relation between R and P. 
It then follows from (29) that R decreases exponentially with the 
time constant 7’. The thermistor voltage U also decreases and for 
small changes this decrease will be proportional to the decrease of 
the resistance. Hence, U also changes exponentially with the time 
constant 7' so that we can obtain 7' from a record of U. (Of course, 
I and U must be so small that the condition N ~ 0 is fulfilled.) 
By determining 7' for a number of different operating points we 
find how 7' depends on P. 


c. Application of the thermistor equation 


An example will be given which shows how the thermistor equation 
(27) can be used for solving a dynamic problem. Let us try to cal- 
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Fig. 11. To the calculation of the transients in the circuit, Fig. 10. 
Nay = + Ra)*s = 


culate U(t) and I(t) in the circuit of Fig. 10, assuming the general 
case where the power N cannot be neglected. 


The heater power is: 
for t | Nir = R,/(R,, + 
for t >0:N,, = 


We assume F#, = constant. Fig. 11 gives the three static U-I- 
characteristics, corresponding to the heater powers N,,, N;2, and 
N, = 0. When the heater power jumps from N,, to N,,, the ope- 
rating point (U, 7) moves along the dynamic characteristic from 
point J to point 2. Our problem is solved if in every moment we know 
the position on the dynamic characteristic of the dynamic point. 
We start at ¢ = 0 in the operating point (U,, /,) and compute the 
operating point at t= 4t < 7. From (27) follows 


At 
AP = (N + — P) (30) 


Here we put N = N, = U, IJ, and P= P, = U;1;, where (Uj, I) 
is the intersection between the resistance line through point J and 
the static U-I-characteristic for N,, = 0. (P has a constant value on 
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a resistance line and can thus be obtained as the thermistor power 
at its intersection with a static characteristic.) We now know 
P=P,+ AP at the time t= At and hence we can draw the corres- 
ponding resistance line in Fig. 10. Its intersection with the dynamic 
characteristic gives the dynamic point at the time t= At. By 
repeating the process a sufficient number of times we obtain the 
solution of our problem. 

Essentially the same method can be used if N;, is a general function 
of time, N;, = N; (t). In eq. (30) we then have to use the average 
value of N,(t) during the time interval in question. 
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CHAPTER 4 


Small Signal Equivalent Circuit 


We shall deduce a small signal equivalent circuit for the indirectly 
heated thermistor. The deduction, which starts from the thermistor 
equation, is based on the fundamental assumption that the variations 
of thermistor temperature and thermistor resistance are so small that 
in performing series developments around an operating point we need 
only consider the linear terms. We further neglect variations in 7' 
and R,, assuming these quantities to have constant values, deter- 
mined by the position of the operating point. , 

We introduce the symbols, given below in (31) and (32). In these 
relations u,7, AR etc. are deviations from the values Uy, J,, Ry etc. 
in the operating point, where the thermistor is in thermal balance. 


N =N, + AN 
I =I, +% (31) Ni, = Nig + 4M, (32) 
R=R,+AR P =P, +47 


The values in the operating point satisfy the relations 
Ry = (33) 
No = (34) 


Introducing (32) and (34) in the thermistor equation (27), we obtain 
d 
AN + 7AN, = AP + To 7 (AP) (35) 


Making use of the power coefficient m, defined in (20), we can put 


AR 
Ry 


(36) 


1 
OP — 
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Here we assume that 


AR < R, (37) 


Introducing (36) in (35) we obtain 
AN AN, : AR+T : AR 
+ = + (AR) (38) 


Evidently, the condition (37) also implies that the change in total 
supplied power, AN + 74AN,, should be small compared to Py. The 
relative changes 4N/N, and AN;/N;,., however, are not necessarily 
small. If, for instance, the heater power increases from N,,. = 0 to 
a value AN, < P,, we obtain an infinite relative change in heater 
power, but eq. (38) is still valid. 

In order to obtain a practically manageable theory, we now assume 
that the variations in current and voltage on the thermistor side 
are small. It follows that we can put 


(39) 


We now introduce (39) into (38) and go over to operational expres- 


d 
sions by putting p = ez’ If we retain the time function symbols 


u, 1, AN, as symbols for the corresponding operators, we obtain 


u i 1 


Introducing in (40) F = — mN, according to (19) and simplifying, 
we find 


1—F-+ pT F nAN,, 


u—iR 0 (41) 
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Fig. 12. Circuit for which equation Fig. 13. Circuit, having the circuit of 
(41) is valid. Fig. 12 as its Thévenin equivalent. 


Evidently, this relation is valid for the circuit of Fig. 12 if in this 


circuit we have 
1—F+ pT 


F = n AN; (p) 


It is easy to see that the circuit of Fig. 12 can be obtained through 
an application of Thévenin’s theorem to the circuit of Fig. 13, if in 
this latter circuit we have 


1—F 
r= (17) 
L=RT 
a) 
en (p) = en (p) (1 + F + pT)/(1 + F) (45) 
Introducing eq. (43) into (45), we obtain 
F 


The power AN, is a quadratic function of current and voltage 
and therefore somewhat inconvenient for our calculations. Provided 
that 4N, < Nj,,5 we can, however, write 4N,, as a linear expression. 
We have N, = J? R, and hence, for small relative changes in N,, 


AN) = 2 = 2 Tron (47) 
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Fig. 14. Equivalent circuit of an indirectly heated thermistor for small superimposed 
currents and voltages. 


Here 7, and wu, are small deviations in heater current and heater 
voltage from the values J,,. and U;,) in the operating point. Introduc- 
tion of (47) in (46) gives us the final expression 


) F Tho 
ep, (p) 1 F U,(P) (48) 


We see that besides the operators e,(p) and u,(p), eq. (48) only 
contains constant quantities, determined by the position of the 
operating point. We can therefore make a direct translation from 
the operators to the corresponding time functions. 

Hence, Fig. 14 gives an equivalent circuit of the indirectly heated 
thermistor for small superimposed voltages and currents in heater 
and thermistor circuit. The quantities r, L, and G are given by the 
following expressions: 


L= RT (49) 


We see that in the equivalent circuit the small heater voltage 
change appears as a small emf in series with the inductance, equal 
to the heater voltage change times a constant G. The value of 
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Fig. 15. Simple circuit studied by the aid of the small signal equivalent circuit of Fig. 14. 


this constant is determined by the position of the operating point 
only. For constant heater conditions we have u, = 0 and come 
back to the equivalent circuit of the directly heated thermistor as 
deduced in Part 1. 

If the heater power in the operating point is so small that the 
condition AN, < Ny» is not fulfilled, we cannot apply eq. (47) or 
use the circuit of Fig. 14. It is easy to see, however, that in this case 
we only have to replace the emf Gu,(t) with an emf G’7AN,(t), 


F 1 


where G’ = 


A pplication 


Let us assume that we have connected an indirectly heated ther- 
mistor to a resistive electric network. By the aid of Thévenin’s 
theorem we can always reduce the network to a series combination 
of an emf FH, and a resistance R,, Fig. 15. Using the equivalent circuit 
of Fig. 16, we can calculate how the thermistor voltage U changes 
at small variations u,(t) of the heater voltage U,. We assume that 
the emf £, is constant. Applying Thévenin’s theorem to the equi- 


i, (t) L r i(t) 


Fig. 16. Equivalent circuit of the arrangement, Fig. 15. ° 
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Fig. 17. The Thévenin equivalent of the circuit, Fig. 16. 


valent circuit, we obtain the circuit of Fig. 17, where the time func- 


tions have been replaced by operators. The impedance Z(p) is 


given by eq. (42). The emf e,(p) follows directly from (43) and (47): 


Further 
€n(P) 
= Zp) + 
and 


u(p) = — i(p) R, 


Making use of (42) we find after some calculations 


Introducing the effective time constant of the thermistor 


t= 
R, — Ro 
we obtain from eqs. (50) to (53) 
2F R, n Ino Up 


(53) 


(54) 


(55) 
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Fig. 18. Another circuit (cf. Fig. 10) for measuring the thermal time constant of an 
indirectly heated thermistor. 


Two cases are of special interest, viz., a) u,(t) changes stepwise, 
b) u,(t) changes sinusoidally. 


a) u,(t) =a sudden step with the height u, 
i From (55) we obtain directly 


| 


i The voltage follows an exponential law with the time constant rt. 
, Kq. (54) gives us 
| 
R= By: t= (57) 
hit 


Thus, when the external resistance is = the thermistor resistance, 
the time constant of the voltage is = the thermal time constant. 

We can use this result for determining experimentally the thermal 
time constant. We employ the arrangement of Fig. 18. Here the 
thermistor forms one branch of a Wheatstone bridge, while the three 
other branches are a, a, and R,. 

First the bridge is brought to balance (R, = R,), so that it has 
the output voltage e, = 0. By short-circuiting the resistance R, < R,, 
we then bring about a small change of the heater voltage. This gives 
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us a small change of the thermistor voltage with a time constant which 
is equal to the thermal time constant (eq. (54) with R, = R, = R,). 
The entire variation of the thermistor voltage appears as the bridge 
output voltage e,(t). Hence, by recording this voltage we can determine 
the thermal time constant. 

A bridge circuit like the one presented in Fig. 18 has two advan- 
tages: we can control in a simple way the condition R, = Ry, and 
we can, by balancing out the constant part U, of the thermistor 
voltage, record accurately the small voltage variation. 

By the aid of the equivalent circuit it is easy to show that we have 
the same time constant for the transient response to a small change 
of the emf on the thermistor side (#, in Fig. 15) as for the transient 
response to a small change of the heater voltage. We can therefore, 
by short-circuiting a small resistance in series with the emf which 
feeds the bridge, also determine 7' from a small change on the ther- 
mistor side. We then come back on a method used in Part 2 for 
determining the thermal time constant of the directly heated ther- 
mistor. In Part 2 the errors of this method are discussed. 

In Chapter 3 we gave another method for the determination of 
the thermal time constant and verified it by a qualitative reasoning, 
based on the thermistor equation. We can now base the verification 
directly on eq. (56). We have assumed that the thermistor power can 
be neglected as compared to the heater power. Hence, F < 1, so that 
(56) can be written 


R, 
wu (t) F R+R 1, u, (1 — e~"*) (58) 
According to (54) we have 
t= (59) 


and thus the thermal time constant is obtained from a record of u(t). 


b) u,(t) = a sinusoidal voltage 


Putting w,(t) = @, sin wt, we obtain the thermistor voltage variation 
u(t) = asin (wt + ~) by replacing in (55) the operator p by jw: 


2F R, n Ino ty 


(60) 


I, 1+ jot 
R, + Ro 
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It follows 


2F Rk, ty 
(61 a) 
R, + Ro 
gy = 1 — arctan wr (61 b) 


We see that the amplitude @ decreases with increasing frequency 
as | / /1 + (w7)?. The value of the phase angle is y ~ 180° at very 
low frequencies and decreases with increasing frequency. At very 
high frequencies it approaches 90°. For wt = 1 the amplitude is 
1 / )2 times its value at very low frequencies while the phase angle 
is 135°. 
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CHAPTER 5 


Experiments 


The measurements described below demanded highly stable ther 
mistor surroundings. In particular, the surrounding medium had to 
be kept at a very constant temperature. The thermistor was there- 
fore enclosed in a metal cylinder, filled with paraffin oil and immersed 
in a temperature controlled bath. 

For some measurements this bath was an oil bath, maintained at 
25° C with a constancy of + 0.01° C by a thermostat with a contact 
thermometer. Particularly for dynamic measurements even such small 
fluctuations turned out to cause unpermissible movements of the 
operating point of the thermistor. Hence, for some measurements the 
thermostat consisted of an ice bath. The metal cylinder with the 
thermistor was immersed in a one litre Dewar flask, filled with 
melting ice made from distilled water. The Dewar flask was placed 
in a cardboard box and thermally insulated with paper. 

The dynamic phenomena were so slow that they could be recorded 
on an xy-recorder (Moseley Autograf, Model 1). The time axis was 
obtained by feeding the x-amplifier with the output voltage from a 
linear potentiometer, driven by a constant speed motor. 

For some of the measurements use was made of an ultralow-frequency 
generator, constructed according to Fig. 19. Its principles are as 
follows. The rotor of a resolver is driven with constant angular speed 
® by a motor working over a gear box. Electrically, the rotor is fed 
from a 1000-cycle A.C. voltage. Through electromagnetic induc- 
tion this voltage gives rise in the stator to a sinusoidally modulated 
1 000-cycle A. C. voltage with a modulation frequency = the rotor 
speed (turns per second). In the resistance R, a small part of the 
stator voltage is added to a small part of the given 1 000-cycle voltage. 
Hence, the input voltage to the amplifier can be written as e’ = 
=A (1+ ksin ot) sin 227-1000t, The amplitude A is regulated with 
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Fig. 19. Ultralow-frequency generator. 


R,, the modulation index k with R,. After amplification the signal 
is demodulated to an output voltage e = HF (1 +k sin wt). A satis- 
factory demodulation requires k < 0.7 to 0.8 and w/22 < 1000 ¢/s. 

A sweep, synchronous with the signal, was obtained by driving 
the time axis potentiometer from the resolver axis with a gear ratio 
1:2. Hence, one period of the signal always occupies one and the 
same length on the record, independently of the «-value. . 


A. Static conditions 


The static characteristics of a number of indirectly heated ther- 
mistors have been accurately determined. These measurements have 
served a twofold purpose: to study the variation of the heater 
efficiency and to investigate the accuracy which can be expected 
from the resistance equation. 

The heater efficiency is of particular interest. Therefore, a meas- 
uring method especially suitable for its determination was chosen. 
In Chapter 2 a general definition of the heater efficiency has been 


given: 
dN 
dN h /}R=constant 


Thus, 7 can be determined from the curves giving N as a function 
of NV, with R as a parameter. 

The circuit of Fig. 20 was employed for obtaining these curves. 
The heater with the resistance 2, forms one branch of a Wheatstone 
bridge; the three resistors a, a, and R, form the other branches. 
From a battery H, the bridge is fed with a current J, which can be 
regulated with the resistor R,. With the switch S in position J the 
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Fig. 20. Arrangement for accurate determination of the static characteristics. 


bridge balance is indicated by the galvanometer G. At balance we 
have R, = R,, giving the heater current J, = + I,. The current J, 
can be accurately determined by the aid of a potentiometer (which 
is not shown in the diagram). The thermistor and the three resistors 
b, b, and R, form a second Wheatstone bridge, similar to the first. 
The second bridge is in balance for R = R, and J = + I,. The current 
I, is also determined in the potentiometer. 

In carrying through a measurement a certain constant value is 
first given to R,. Then, with J, = I, = 0, the value of the current I 
which gives R = R, is determined (bridge 2 in balance). J is then 
decreased and both bridges brought to balance simultaneously, bridge 
1 by regulating R,, and bridge 2 by varying the heater current. 
At the same time the two currents J and J, are also determined. 

A pair of corresponding values of the thermistor power N = RI? 
and the heater power N, = R, 1; has now been obtained. The 
current J is further diminished and a new pair of values (N, N,) is 
determined and so on. The measurements also display the variation 
of R,. 

An investigation of several different thermistor types gave entirely 
consistent values. Fig. 21 shows the results for the thermistor of 
Figs. 4 to 9 (Stantel, type B 2552/60) at 0, = 0°C. The diagram 
gives N/P vs. N,/P, where the dissipated heat power P has been 
determined as the thermistor power corresponding to N, = 0. Evi- 
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(9) 0.2 0.4 0.6 0.8 1.0 1.2 


Point 1-2 3-7 8-9 0-11 12-15 16-19 
Rk 100 40 10 4 1 0.4 
PoW 5.97 9.96 17.39 23.57 35.96 47.23 


Fig. 21. Relation between N, Nj, and P for the thermistor of Fig. 4. 


dently, all measured points are very close to a straight line. Hence, 
we can write with good accuracy the relation between N, N,, and R as 


N +N, = P(R) 


with » = constant and P(R) = the thermistor power which for 
N;, = 0 gives the thermistor resistance R. From Fig. 21 we find 
n = 1/1.135 = 0.881. 

Only two of the measured points, 18 and 19, deviate noticeably 
from the straight line. These points correspond to relatively high 
thermistor temperatures. All the other measured points show a devia- 
tion from the straight line which corresponds to a variation in 7 
of only a few tenths of a percent. This small variation is entirely 
covered by the measuring errors. 

Measurements at 0, = 25°C gave the same value of 7 as those 
at 0,=0°C. This indicates that 4 is also independent of the 
ambient temperature. ; 
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Fig. 22. Relation between R, and P for the thermistor of Fig. 4. A straight line 
approximation is suggested. 


In calculating N, = R, 1% one must pay regard to the variation 
from point to point of R,, as illustrated by Fig. 22. The measured 
points are all within the hatched part. The lower border line corre- 
sponds to NV, = 0, the upper border line to N = 0. Hence, strictly 
speaking, two different pairs of values (N, N,), which correspond to 
one and the same value of P, give somewhat different values of 
R,. In other words, the relation between R, and P = N + 7N, 
is not an exact one-to-one relation. 

In fact, if we increase the heater power and, at the same time, 
decrease the thermistor power so as to keep P constant, we obtain 
a small increase of R,. The reason is that at large heater power and 
small thermistor power the entire heater wire becomes heated, while 
at small heater power and large thermistor power only that part 
of the heater wire becomes heated which is in direct contact with the 
thermistor body. In most cases this variation can be neglected so 
that the relation between R, and P can be approximated by a straight 
line as indicated in Fig. 22. 

The resistance equation (10) applies with excellent accuracy to 
the thermistors investigated. Consider, for instance, the thermistor 
B 2552/60, mentioned above. The table in Fig. 21 gives six points 
(R, P) on the R-P-characteristic. From three of them (100 kohms; 
5.973 mW), (10 kohms; 17.39 mW), and (1 kohm; 35.96 mW) 
we compute the three constants R,, B, and C. This gives us the 
resistance equation 
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From this expression the resistances in the three other points are 
derived with a maximum error of about 1 percent. The same 
accuracy is obtained if we use the equation for calculating points on 
the R-P-characteristic for 0, = 25° C. 


B. Dynamic conditions 


The dynamic measurements were performed on a thermistor, type 
B 2321/60, which was kept at the temperature 0° C (OQ) = 273° K). 
The resistance equation at this ambient temperature came out as 


2960 
R = 0.0931 e558? + 273 ohms 


The heater efficiency was 7 = 0.862 so that 
P=WN 0.862 


Figs. 23 and 24 show the R-P-characteristic and the static 
U-I-characteristics. 
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Fig. 23. R-P-characteristic of the thermistor Stantel type B 2321/60 at 0°C. 
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Fig. 24. Static U-I-characteristics of the thermistor of Fig. 23. 


a) Nz = 0 d) N, = 30 mW 
b) N, = 10 mW e) N;, = 40 mW 
c) N, = 20 mW 


Eq. (21) gives the dynamic factor: 


530 N 
~ (P+ 48.9) 


F has its maximum value F,,,, = 2.71 for N,=0, N= P= 
= 48.9 mW. 

The heater resistance varies between 98.8 ohms (for P = 0) and 
101.2 ohms (for P = 50mW). The average value R, = 100 ohms 
has been used. 


a) Stepwise variations, thermal time constant 


The thermal time constant has been determined by the aid of the 
method of Chapter 4. In principle the bridge circuit of Fig. 18 was 
used, where transients can also be excited by short-circuiting a small 


x 
T\ 
i 
| 
| 
| 
i 
! 
i 
| 


39 
T sec 
12 

8 
4 
0 P oW 

ie) 10 20 30 40 50 


Fig. 25. Thermal time constant 7 as a function of P for the thermistor of Fig. 23. 


resistance inserted in series with H,. From about 30 measurements at 
different operating points the following conclusions have been drawn: 


1) Changes on the thermistor side and on the heater side give the 
same time constant. 

2) The thermal time constant is practically the same for different 
values of N and N, as long as P = N + 7N, is kept constant. 
There is possibly a tendency that 7' increases with increasing N,,. 
The observed variations, about 5 percent, were of the same order 
of magnitude as the measuring errors. 


The curve of Fig. 25 gives the relation between 7 and P yielded 
by the measurements. In all measured points the deviation from the 
value given by this curve was less than 5 percent. 

Some determinations of 7’ at very small thermistor powers were 
made according to the method described in Chapter 3. The results 
were in complete agreement with those obtained by the bridge method. 
The circuit of Fig. 26 was employed. The emf HL, compensates for 
the thermistor voltage in the operating point so that only the voltage 
variations are recorded. The emf £#, and the resistance R, were 
chosen so as to give N ~ 0 and thus F < 1. 

According to eq. (54) the time constant of the recorded voltage 
becomes equal to the thermal time constant if F be <i. 

R, + Ro 
This condition is evidently fulfilled if R, ~ R,. As we have chosen 
F < 1, we can, however, tolerate a fairly large difference between 
R, and R, without jeopardizing the accuracy. This is an advantage 
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Fig. 26. Arrangement for recording small variations of the thermistor voltage. 


against the bridge method where already a small difference between 
R, and R, gives a large error, so that it becomes necessary to correct 
for the influence of either the input impedance of the recording 
instrument or the resistance of the battery branch (cf. Part 2, p. 23). 
In most cases, therefore, the method of Fig. 26 is to prefer. 

Table 1 gives the result of measurements for such values of the 
series resistance R, (see Fig. 15) that the effective time constant 1 of 
the thermistor differs considerably from the thermal time constant. 
The circuit of Fig. 26 was employed. Considering the many possible 
sources of error, the agreement with theory must be considered as 
satisfactory. 


TABLE 1. Effective time constant t. Comparison between theoretical values from 
eq. (54) and experimental values, obtained from transients. 


R, P N Ry T sec. 
ohms mW sec. | mW mW ohms theory | exper. 
500 13.0 ll 0 13.0 1.80 2 500 5.0 5.2 
10 4.4 0.61 2 500 7.9 i 
10 4.4 0.61 100 18 15.6 
50 35.6 10 20 18.4 1.37 250 5.3 5.6 
40 1.2 0.09 4 000 9.1 wa, | 


b) Sinusoidal variations ; 

With the aid of the ultralow-frequency generator, described above, 
an experimental determination was made in different operating points: 
1) of the thermistor impedance for small superimposed sinusoidal 
voltages on the thermistor side, 2) of the frequency dependence of 
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Fig. 27. Determination of the thermistor impedance for small superimposed sinusoidal 
voltages. Fig. 19 shows the principle of the ultralow-frequency generator. 


a small variation of the thermistor voltage, caused by a small sinusoidal 
voltage superimposed on the heater voltage. 

The output voltage of the ultralow-frequency generator was of the 
form e = H+ k Esin wt, where E and k could be varied indepen- 
dently of each other. At first, with k = 0, H was chosen so that 
the desired operating point was obtained. Then a small A.C. voltage 
was superimposed by regulating & to a suitable value, generally 5 to 
10 percent. 


The circuit of Fig. 27 was employed for the impedance measure- 
ments. The y-amplifier of the xy-recorder was supplied with the 
thermistor voltage variations. The thermistor current variations 
caused variations of the voltage over R, and these variations were 
brought to the x-amplifier of the recorder. The two variable emf’s 
connected in series with the input to the recorder compensate for the 
constant parts of voltage and current. For each frequency the rela- 
tion between voltage and current was recorded. From the ellipse 
obtained the absolute value and the phase angle of Z were easily 
calculated. 

According to eq. (42) the thermistor impedance can be written 


1—F+joT 
°14+F4+j0T 


Z(w)=R 


In the complex plane Z is represented by a half circle (see Parts 
1 and 2). For o=0 we have Z= R, (1 — F)/(1+ F) =r; for 
«w = » we have Z = R,. In the range F > 1 the phase angle of Z 
goes from 180° to 0° as w goes from 0 to. 

Fig. 28 gives theoretical and experimental impedance diagrams for 
a number of operating points, corresponding to one and the same 
value of the static resistance. As in Part 2, for the directly heated 
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Fig. 28. Impedance diagram for the thermistor of Fig. 23. In all cases the static 
resistance of the thermistor is 50 ohms. 1 me/s = 0.001 ¢e/s. 
a) Ny =.0 ; N = 35.6 mW 
b) Nz = 10 mW; N = 27.0 mW 
c) Nz = 20mW; N = 18.4 mW 
d) = 30mW; N= 9.8 mW 


thermistor, we find that an experimentally determined curve always 
lies inside the corresponding theoretical curve. The experimental 
value of the phase angle is always smaller than the theoretical 
value. Mostly the absolute value of the impedance is also smaller 
than the theoretical value. At low frequencies it can be larger. 


The arrangement of Fig. 29 was used for recording the small 
thermistor voltage, caused by the superposition of a small sinusoidal 
voltage on the heater side. At each frequency two cycles of the 
superimposed thermistor voltage were recorded, and from the records 
the amplitude and phase characteristics were determined. 

According to theory (cf. eq. (60)) the superimposed thermistor vol- 
tage should vary as 1/(1-+ jr). This means that: 1) at the critical 
frequency /,, determined by 2z7f,t = 1, the amplitude has decreased 
3 dB, the phase angle 45°, 2) at high frequencies the amplitude varies 
as 1/wr, i.e., it decreases 6 dB per octave. 
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Fig. 29. Arrangement for determining the amplitude and phase characteristics when 
superimposing a small A.C. voltage on the heater voltage. 


Fig. 30 shows three pairs of amplitude and phase characteristics, 
corresponding to three cases with the same operating point but 
different values of the series resistance R,. The curves are the theo- 
retical characteristics, corresponding to the values of f, stated in the 
figure. The measured points are very close to the curves. At all 
measurements performed this was found to be the case. 
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Fig. 30. Amplitude and phase characteristics for the thermistor of Fig. 23. 
Same operating point in all cases. 
1 me/s = 0.001 e/s. Ry = 500 ohms. 
a) Ry = 100 ohms; b) Ry, = 500 0hms; c) Ry = 2 500 ohms. 
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Table 2 gives the results of three different series of measurements. 
The experimental value of t has been determined from the critical 
frequency, t = 1/2zf . 

In Table 2, Series No. 1 covers the measurements described above, 
performed in one and the same operating point but with different 
values of R,. The agreement between theory and experiment is 
good. Series Nos. 2 and 3 correspond to different operating points 
with the same static resistance but with different heater powers. In 
Series No. 2 the series resistance R, is equal to the static resistance, 
in Series No. 3 the series resistance is much larger. For R, = Ry, 
the agreement is very good in all cases, for R, > R, there are 
appreciable deviations in those cases where the thermistor power is 
large. 

A closer study of the experimental results of Tables 1 and 2, and 
of the experimental impedance diagrams, shows that the deviations 
between theory and experiment can be explained by the hypothesis 
that the value of the dynamic factor F is different under static and 
dynamic conditions. In other words: the dynamic resistance of the 
thermistor is different in static and dynamic (time-varying) pheno- 
mena. 

If the experimental results are used for calculating the “correct” 
value of F, we find on the one hand that this “correct” value is 
always smaller than the theoretical value, on the other hand that 


TABLE 2. Effective time constant t. Comparison between theoretical values from 
eq. (54) and experimental values, obtained from the sinusoidal state. 


| Ry P Ni, N Rp T sec. 
| ohms mW sec mW mW ohms theory | exper. 
500 13.0 11 10 4.4 0.61 100 18 16,8 
500 11 11.2 
2 500 7.9 9.0 
50 35.6 | 10 10 27.0 2.01 50 10 10.0 
20 18.4 1.37 10 10.1 
30 9.8 0.73 10 10.0 
40 1.2 0.09 , 10 10.1 
50 35.6 10 10 27.0 2.01 4 000 3.4 4.5 
| 20 18.4 1.37 4.3 5.2 
30 9.8 | 0.73 5.9 6.8 
40 1.2 0.09 9.1 9.2 
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it depends both on the frequency and on the series resistance R,. 
If we turn to the original definition of the dynamic factor, 
F Gea it foll that rtain ch in th 

GN cons.’ it follows that a certain change in the 
thermistor power during a dynamic course of events always gives a 
smaller resistance change than expected. Of course, this is due to 
the fact that we have neglected that under dynamic conditions the 
temperature distribution is uneven. 

It is evident that this uneven temperature distribution must depend 
both on the frequency and on the properties of the circuit in which 
the thermistor is inserted. The measurements agree well with theory 

R, +R, is < 1, i.e., where R, 
is about equal to R, or where the thermistor power N is so small 
that F ~ 0. 


in those cases where the product F 
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